Positrons in Cosmic Rays from Dark Matter Annihilations 
for Uplifted Higgs Regions in MSSM 
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We point out that there are regions in the MSSM parameter space which successfully provide 
a dark matter (DM) annihilation explanation for observed e + excess (e.g. PAMELA), while still 
remaining in agreement with all other data sets. Such regions (e.g. the uplifted Higgs region) 
can realize an enhanced neutralino DM annihilation dominantly into leptons via a Breit-Wigner 
resonance through the CP-odd Higgs channel. Such regions can give the proper thermal relic DM 
abundance, and the DM annihilation products are compatible with current antiproton and gamma 
ray observations. This scenario can succeed without introducing any additional degrees of freedom 
beyond those already in the MSSM. 
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INTRODUCTION 



There has been growing interest in the interpretation 
of an observed positron excess in light of the recent data 
from PAMELA 0-Q and Fermi-LAT 0-11], while satis- 
fying antiproton and gamma ray constraints. 

Among the possible sources for the positron excess are 
astrophy_sical sources, such as pulsa rs [ 7H12| and super- 
novae [iMLzl], and dark matter fT8j-[43|7^ The dark mat- 
ter (DM) possibility is of great interest from the particle 
theory viewpoint, and we seek in this paper the cosmic 
and gamma ray signatures of DM annihilation in a su- 
persymmetric model within the MSSM (Minimal Super- 
symmetric Standard Model) parameter space. 

Most of the conventionally explored MSSM parameter 
space cannot successfully explain the observed positron 
data. There is typically a significant branching ratio of 
DM annihilation into gauge bosons and Higgs as well 
as quarks and hence too much hadronic production of 
antiprotons, in excess of what is observed. The MSSM, 
however, has more than a hundred parameters and it 
would be worth seeking such a possibility to realize the 
current cosmic ray observations in the framework of the 
MSSM without introducing any additional degrees of 
freedom. Another persistent problem in the dark matter 
annihilation scenarios to explain the observed positron 
excess is the requirement of large annihilation cross sec- 
tions in the halo far bigger (typically a boost factor of a 
factor 100 or more) than the canonical thermally aver- 
aged cross section value at freeze-out for the weakscale 
dark matter inferred from the observed dark matter relic 
density. Several resolutions to this apparent discrep- 
ancy in the annihilation cross sections have been pro- 
posed, such as unconventional cosmological histories |4# - 
|46| which can affect the dark matter freeze-out temper- 



ature; Sommerfeld enhancement [47H49l | which requires 
new light particles to allow new long range interactions 
between dark matter particles; Breit-Wigner enhance- 
ment [50l - l53j which requires a particle whose mass is 
close to the twice of the dark matter mass; and sub- 
structure clumps (54l - [59| which could provide a partial 
contribution to the enhancement (say, by a factor of a 
few). The annihilation enhancements make the the cur- 
rent annihilation cross section higher in galaxies today 
than it was at the time of freeze-out, so that it can ex- 
plain the positron excess (which requires a high cross 
section today) while still satisfying the correct relic den- 
sity (which requires a lower cross section at freeze-out). 

The main results of the paper are to point out the 
existence of parameter regions in the MSSM which can 
potentially realize a DM annihilation scenario that ex- 
plains the observed positron excess. We show that there 
can be dominant leptonic final states and a large boost 
factor to obtain the positron excess and thermal relic 
abundance without introducing any additional degrees 
of freedom beyond those already in the MSSM; in addi- 
tion these scenarios can be compatible with the current 
antiproton and gamma ray data. 

As a concrete example, we consider the following sce- 
nario: The requirement of a dominant leptonic final state 
(in this case taus) can be satisfied in the uplifted Higgs 
regions [60j within the MSSM as described in Section [ft] 
In addition, we obtain the required boost factor from the 
pseudo-scalar Higgs s-channel resonance in the MSSM 
which can induce the Breit-Wigner enhancement [50l — 
53| via XX ~ * A — > t + t~ for tua ~ 2m x (x denotes 
the neutralino dark matter and A denotes the CP-odd 
Higgs). 

We also note that we have here found a situation where 
annihilation to taus satisfies all existing observations. 
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The reason is that we treat the FERMI data as astro- 
physical background (fit by a simple power law) while 
we take the PAMELA excess to be due to DM annihi- 
lation. Our interpretation is in contrast to some of the 
previous literature, where the requirement was made of 
explaining the excesses in both PAMELA and FERMI 
as due to DM annihilation. 

In this paper we take the dark matter to have a cored 
isothermal density profile [HJ, |(32| in the Galactic Halo 
unless stated otherwise. This relatively flat core density 
distribution helps to ameliorate the severe gamma ray 
constraints |63l465| (from observations e.g. by FERMI 
and HESS) in contrast to other profiles such as the NFW 
profile 66!] . Our choice of profile should be adequate for 
the purpose of illustrating the potential significance of 
the previously unexplored MSSM parameter space dis- 
cussed here. 

Sec. HTl reviews the uplifted supersymmetric Higgs re- 
gions in the MSSM parameter space and its unique prop- 
erties well motivated for the current cosmic ray observa- 
tions. Sec. Mil then discusses the cosmic and gamma 
rays signals expected for such regions, followed by the 
conclusion/discussion in Sec. [IVJ 



II. MOTIVATION: UPLIFTED HIGGS REGION 

The purpose of this section is to briefly review the 
uplifted Higgs regions (heavy Higgs mass is 'uplifted' in 
those regions) and point out that there exists a region 
in the MSSM parameter space where the dark matter 
can annihilate dominantly to leptons and its annihila- 
tion cross section can be boosted via the Breit-Wigner 
resonance without introducing any additional fields be- 
yond those already in the MSSM. We refer the readers 
to Ref. [6(| for more detailed discussions on the uplifted 
Higgs regions. 

In the usual MSSM, it is difficult to obtain the re- 
quired leptophilic cross sections to explain PAMELA 
while not overproducing p. In particular, in this pa- 
per we focus on the case where the dominant channel 
for DM annihilation takes place via Higgs resonance 
XX A -> T + r~,bb. In the standard MSSM, for those 
Higgs funnel regions with a large tan /?, the problem is 
that the branching ratio to t+t~ is typically at the 10% 
level while the remainder is predominantly to bb; the lat- 
ter produces far too many p in excess of what is observed. 

In the standard case the ratio of final state t + t~ to 
final state bb must be small because BR(xx — ► A —> 
t + t-)/BR(xx -> A ->• bb) ~ yl/Zyl and the Yukawa 
couplings y b > y T since m&(= yb(H d )) > m T (= y T (H d )). 
In the uplifted Higgs regions, on the other hand, the 
bottom type fermion mass and the bottom type Yukawa 



coupling are not necessarily proportional to each other, 
so that one can have y T > yb and the dominant an- 
nihilation to t. Unlike the standard case, here m T is 
generated not by the usual tree level down-type Higgs 
vacuum expectation value (VEV) (Hd) (which vanishes 
at tree level) but instead at the one loop level from (H u ) 
and (Ha). 

Such unique features in the uplifted Higgs regions 
arise by the absence of the soft-SUSY breaking B-term, 
BjiH u H c i (which can be justified for instance by the 
appropriate R-charge assignment R(Hd, Q, U c , E c ) = 
and R(H U , D c , L) — 2) while the superpotential is same 
as that of the conventional R-parity conserving MSSM. 
This results in the vanishing down-type Higgs VEV at 
the tree level by enforcing the electroweak symmetry 
breaking via \n\ 2 + m 2 Hu < 0, \n\ 2 + m 2 Hd > and the sta- 
bility of H u Hd flat direction via 2|/i| 2 + 771^ +m 2 Hd > 0. 
(Hd) ^ comes from the loop contributions and it helps 
in explaining the small down type masses which are loop 
suppressed. The down type quark and lepton masses 
arise at the one-loop level 

m d = y d (H d ) +y' d (H u ) (1) 

yd is the standard down type Yukawa coupling between 
Hd and down type fcrmions and y' d represents the loop- 
induced effective Yukawa coupling between H u and down 
type fermions [60S]. Hence yd is not in general pro- 
portional to nid in the uplifted Higgs regions anymore, 
and, for instance, y T > yb can be possible even though 
m r < nib- 

In particular the decay of a heavy Higgs into the 
down type fermion pairs is proportional to yd but not 
to nid tan j3 in contrast to the usual MSSM. The param- 
eter space we focus on has tua ~ 2m x and tan (3 ^> 1 as 
well as y T > yb, so that, as discussed above, we consider 
the dominant dark matter annihilations through Higgs 
resonance XX ~^ A — > t + t~, bb (note A is dominated by 
Hd component in our scenarios). The coupling of A to 
the top quarks are highly suppressed by the factor cot 2 /3 
and the final states are dominated by r and b with an 
aforementioned ratio Br(xx T+T ~)/Br(xx ~* bb) ~ 
y 2 /(3y 2 ) in our scenarios. Such regions have been missed 
in the conventional studies of MSSM, and we argue the 
potential significance of such regions for its unique cos- 
mic ray signals distinguishable from the conventional 
MSSM cosmic ray predictions. 

We in the following briefly discuss the enhancement 
of the annihilation cross sections in the uplifted Higgs 
regions to see the range of the parameters required for a 
sufficient boost factor . We simply assume the dark mat- 
ter is a neutralino LSP (lightest supersymmetric parti- 
cle) and we are interested in the case where its mass is 
close to the twice of the pseudo-scalar Higgs to obtain a 
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large boost factor via the Higgs resonance. The relevant 
couplings in the Lagrangian are 



with 



C 



L 3 C ddA d lb dA + ±C™ A X75XA 



V2 



(2) 



C xxA = -i(gN 12 -g'N n ){smf3N 13 -cospN u )(3) 

where Nij is the neutralino mixing matrix in the 
(B,W°,H%,H°) basis, g',g are U(1),SU(2) gauge cou- 
plings respectively and our scenarios are dominated by 
d = r,b final states. The dark matter annihilation cross 
section through the Higgs resonance a xx (xx ~ ► A —> dd) 
becomes 
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where s is the Mandelstam variable and p — 



s — 4m|/2 represents the initial momentum of the in- 
coming particle x m the center-of-mass (cm.) frame, 
k = yfs — Am 2 /2 is the final cm. momentum and 
is the total decay width of A. The branching ra- 
tio Br(xx ^_bb) ~ O(0. 1) with Br( X x -> t + t~) + 
BR(xX ~ * bb) « 1 can be possible in the uplifted 
Higgs regions in contrast to the usual Higgs resonance 
in MSSM with a large tan/3 where Br(xx t + t~) ~ 
O(O.l). Its thermal average is [671168} 
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6), 7 = T A /m A - 
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(7) 



The conventional calculation of the neutralino annihila- 
tions in terms of the series expansions of x is inadequate 
when m,A ~ 2m x and Figure [1] shows the numerical inte- 
gration of {av) normalized by its value at the freeze-out 
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FIG. 1: Boosted cross section: The numerically in- 
tegrated Breit-Wigner enhancement {av) (x) / {av) (x/ = 
25) for, from the top to bottom curves, <5 = 7 = 
10- 5 ,10- 4 ,10- 3 ,10- 2 ,10- 1 as denned in Eq.Q. {av){x = 
25) represents the usual thermally averaged annihilation cross 
section at 'freeze out' x = m x /T = 25. 



temperature Xf (we set Xf = 25 for illustration) for a 
few parameter sets (S, 7) . 

Compared to the non-resonant case with constant (av) 
after the freeze-out, the dark matter abundance with 
the Breit-Wigner resonance is enhanced by the effective 
boost factor (BF) ^ 



BF 



-*oo, no 



{<TV)\T=0/Xf 

f°°dx^ 

JXf X z 



(8) 



where we took the same cross section at zero temperature 
to be (<jv)\t=o- The present relic abundance hence is 

^,Q.lx 1Q ; 9 f r2 ^x^ (9) 
x M| T =o v 7 ^ 

where is the effective number of degrees of freedom. 

We compute values for the boost factor using 
Eqn.([8|). For instance, the choice of S = 7 = 1CP 1 
(10~ 2 ,10~ 3 ,10~ 4 ,10~ 5 ) gives respectively the boost fac- 
tor of the order a few (0(10), 0(50), 0(300), 0(3000)). 
The justification of these small S, 7 values is beyond the 
scope of this paper. Nonetheless, we consider the exis- 
tence of such leptophilic boosted regions to be a strong 
enough motivation for us to seek compatibility of this 
MSSM parameter space with recent cosmic ray observa- 
tions. 

In the following sections, we focus primarily on the 
dark matter mass of the electroweak scale. We allow the 
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branching fraction Br(xX t + t~) to vary up to 90% 
with Br(xx — > t + t~) + BR(xx — * bb) « 1 as expected 
for the Higgs resonance regimes in the uplifted Higgs 
regions [60|. This study should suffice for the purpose 
of pointing out the potential significance of previously 
missed MSSM parameter regions in view of the recent 
cosmic ray data. 

III. COSMIC RAY SIGNALS 

In this section we discuss positron, antiproton, and 
gamma ray signals in our scenario. Henceforth for 
the density profile of the gala ctic halo we take the 
cored isothermal profile [6l|, [62| to alleviate the strin- 
gent gamma ray constraints as we shall see, and we use 
the MED model [jO, [7l| for the propagation parameters 
unless stated otherwise. We also use the force field ap- 
proximation [72l [73| with the solar modulation parame- 
ter (Fisk potential) 4>f = 0.6 GV for definiteness for the 
charge independent solar modulation effects and do not 
consider the charge dependent solar modulation [l8|, 
whose realistic treatment is left for future work due to 
large uncertainties. 

A. Positron signals 

We attempt to explain PAMELA e+/(e+ + e~) data 
as the signature of dark matter annihilation. However, 
we treat FERMI e + + e~ data as background; then the 
DM additional contribution to e + + e~ must be at most 
a small addition on top of this background. 

Positrons in the Galactic halo travel under the influ- 
ence of the interstellar magnetic fields and lose energy 
via inverse Compton scattering and synchrotron radia- 
tion. Such effects (essentially random walk processes) 
can be treated via the steady-state diffusion equation 
for the number density per unit energy N(E, x) |75l - l78| 

^ = = V ■ [K(E, x) ■ V A] + A [ h {E, k)N] + Q{E, x) 

where K is the diffusion coefficient and b = (1/r) x 
E 2 /(lGeV) is the energy loss rate (we take r = 10 16 
seconds) and Q is the source term which is, in the dark 
matter annihilation scenarios, proportional to the dark 
matter number density squared and the thermally aver- 
aged annihilation cross section. 

First we discuss our treatment of the FERMI e + + e~ 
data as background. A simple featureless power law, for 
instance the one in Ref. [33 
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can fit the FERMI e + +e~ data as shown in Figure[2](thc 
statistical and systematic errors are added in quadra- 
ture). We are in this section primarily interested in 
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highest curve, which is the sum of the DM contribution plus 
the FERMI fit, is still within the systematic errors of the 
FERMI data. 

the dark matter mass of order a few hundred GeV and 
the boost factor ~ 0(100) which can affect the lower 
energy part of FERMI data by order ~ 0(10)% as il- 
lustarted in Figure [2] Such a 10% shift in the spectrum 
can well be within the systematic errors [37|, and we do 
not expect, in the scenarios of our interest in this paper, 
the dark matter annihilations to affect significantly the 
current FERMI e + + e~ spectrum. 

Now we turn to PAMELA data. In contrast to the 
relatively smooth FERMI e + + e~ spectra, the bump- 
like rising behavior of the positron excess by PAMELA 
would be hard to represent merely with a simple power 
law, and we interpret this anomalous positron excess as 
the signal of dark matter annihilations. 

For the astrophysical background of the positron frac- 
tion, we make a simple featureless power law fit for the 
PAMELA low energy data set 



0.07 
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which was conservatively obtained from the lowest 6 en- 
ergy bins out of the total 16 bins in the PAMELA data 
recently analyzed in Ref. without taking the charge- 
dependent solar modulation (whose uncertainty is quite 
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FIG. 3: The uncertainties in the positron fraction back- 
grounds. The most recent PAMELA data have a some- 
what lower positron fraction than the earlier analysis. As 
our backgrounds, we use the second curve from the bottom 
(the pink dotted line), which is an extrapolation of the low- 
energy PAMELA data and lies in between the widely-used 
higher Moskalenko and Strong (MS98) result and a fit with 
the spectral index 7 = 3.44 (labeled med electron). 
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FIG. 4: The electron flux background $ e - . 



large) . We extrapolate this fit to higher energies and take 
it to be the e + /(e + + e~) background. This fitting for- 
mula gives a background higher than the one convention- 
ally used (as we discuss shortly). Including the higher 
energy bin data in order to get a fit makes the back- 
ground even bigger. A few comments on this positron 
fraction background are in order. Figure [3] illustrates 
the uncertainties in the positron fraction background es- 



FIG. 5: The positron fraction e + /(e + + e ). The branch- 
ing ratio of DM annihilation to taus is as labeled, where 
BR-(XX -» t + t~)+BR(xx ->■ bb) = l. Boost Factors for the 
different curves are as labeled. 



timations. Our background obtained by a simple ex- 
trapolation of PAMELA data is bigger than the conven- 
tionally used positron fraction background obtained by 
the electron and positron fluxes of Moskalenko & Strong 
(MS98) [3, [[3. The widely used MS98 fits can po- 
tentially underestimate the positron fraction partly be- 
cause of the overestimation of the electron flux [80j as 
shown in Figure 2J which contrasts the MS98 fits with 
the more recent electron flux data points from CAPRICE 
and AMS [U |82j]. Figure 0] also shows the electron 
flux parametrized by a power law with spectral index 
7 = 3.44 ± 0.1 [HI, which better represents the re- 
cent electron flux data and gives a smaller e~ amplitude 
than that used by Moskalenko & Strong. The three in- 
dices 7 = 3.34, 3.44, and 3.54 are labeled "hard elec- 
tron", "medium electron" and "soft electron" respec- 
tively in the figure. The positron fraction background 
obtained by combining the electron flux with 7 — 3.44 
(labeled med) together with the positron flux estimation 
by Moskalenko & Strong are also shown in Figure [3] The 
current data of the positron flux itself still have too large 
error bars for parametrization fitting to be feasible; how- 
ever the forthcoming positron flux data such as those 
from AMS-02 84] would clarify these issues. Figure |3] 
also shows the previous PAMELA analysis [l[ (denoted 
as old PAMELA) for reference. For definiteness, in the 
following, we use the most recent PAMELA analysis @ 
(denoted simply as PAMELA) and the simple PAMELA 
low energy data extrapolation for our positron fraction 
background. Our choice of background seems reason- 
able because it lies in between the conventional MS98 
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positron fraction background and the alternative back- 
ground obtained by using MS98 e + together with the 
"medium electron" fit to the currently available e~ data 
as can be seen in Fig [3J 

Figure [5] shows the PAMELA positron fraction ex- 
cess data together with the predicted positron fraction 
from DM annihilation in our scenarios. One can see 
that the following parameter choice is a good fit to data: 
the boost factor BF=110, BR(xx -> r + r - )=0.9, and 
BR(xX —> 66)=0.1. For this parameter choice, the av- 
eraged x 2 for 10 highest energy bins is 1.0 and that for 
the 8 highest energy bins (the choice of these bins is as 
conventionally chosen in the literature to avoid the ef- 
fects of the charge dependent solar modulation) is 0.7. 
Since this set of parameters is a good fit to PAMELA 
data, we take it to be our canonical case for upcom- 
ing discussions of other cosmic ray signals. Indeed this 
choice will illustrate the existence of the MSSM param- 
eter space compatible with the current cosmic ray data. 
We note that, had we chosen a different background fit 
with a bigger amplitude such as the one shown in Figure 
|3J a good match of our predictions to the PAMELA data 
would be even easier because the required contribution 
from DM would be smaller (e.g. the boost factor could 
be smaller). 

In contrast, for the case with DM annihilation primar- 
ily to 66, the match to PAMELA data is quite poor. This 
is the typical situation in the usual MSSM. In the figure 
we illustrate the case with BR(xx — ^ t + t~ )=0.1 and 
BR(xX ~~ ► 66)=0.9; here the curves shift toward lower 
energy because of the softer positrons from 6 compared 
with those from tau decays. 

One can obtain a harder positron fraction spectrum 
by choosing a higher dark matter mass (the cutoff scale 
of the positron fraction spectrum shifts towards higher 
energy for kinematic reasons) and an example for m x — 
400GeV is shown for comparison. The higher dark mat- 
ter mass decreases the overall amplitude due to the 
smaller number density, and hence we compensated for 
it with a bigger boost factor 200 leading to the averaged 
X 2 = 1.0 for 10 highest energy bins and x 2 = 1A for the 
8 highest energy bins. 

Even though the exclusive parameter scanning for the 
uplifted Higgs regions in the MSSM parameter space is 
beyond the scope of this paper due to the nontrivial par- 
ticle physics constraints (such as those from flavor chang- 
ing interactions), we can see that a wide range of MSSM 
parameters in the uplifted Higgs region can be consistent 
with the current positron fraction data. 

Now we turn to constraints from other cosmic ray data 
for consistency. 



B. Antiproton Flux 

The antiproton flux could be problematic because 
our scenario expects the 66 annihilation channel to be 
typically of order 10% or more, and the non-negligible 
hadronic processes can overproduce antiprotons in excess 
of what is observed. The antiproton flux is obtained by 
solving the steady diffusion equation [7(1 HE| analogously 
to that for electron/positrons (except for the negligible 
energy loss term because m p m e ) 

^— = = K{T)V 2 N - ^-(stgn(z)NV c ) (11) 
at oz 

+Q{*,T)-2h8{z)Y ann N 

where T is the antiproton kinetic energy, V c is a constant 
galactic convective wind (corresponding to a constant 
flow of magnetic irregularities) and h = O.lkpc is the 
height of our Galaxy approximated as a thin disk, T ann 
denotes the annihilation rate of antiproton and the inter- 
stellar proton. Even though there still exist large uncer- 
tainties in the background estimations of the antiproton- 
to-proton ratio, we present the following discussions us- 
ing the secondary antiproton background estimated by 
the fitting formula [86] 

log 10 4>* G = 0.028log 4 w (T/GeV) - 0.02log 3 w (T/GeV) 
-1.0 log 2 (T/GW) +0.071og 10 (T/GeV0 - 1.64 

and that for the primary protons via [87| 

j.BG/rp\ = 890(r/GW)-"- 43 

9p 1 ' 1. - 0.112(T/GeVO°- 992 + 0.156(T/GeF) 2 - 03 

Figure [6] shows the dark matter contributions to the 
antiproton-to-proton ratio for m x = 300 and BF=110, 
where we conservatively treated the PAMELA data Q 
as the upper bound because of the large uncertainties 
in the background estimations. Figure [5] clearly illus- 
trates the advantage of the uplifted Higgs regions where 
the branching fractions into taus can be quite large com- 
pared with other conventionally studied MSSM param- 
eter regions with a large tan ft. In fact, the conventional 
MSSM parameter space is typically not compatible with 
the lack of antiproton excess in the PAMELA. Indeed the 
cases with branching ratio into 66 well over a few tens of 
percent are strongly constrained. Taking advantage of a 
small branching fraction to 66 of > (9(10)% in our dark 
matter annihilation scenarios, we can see that our sce- 
narios can be compatible with the lack of the antiproton 
excess. 

For illustration purposes, we also added the back- 
ground contributions to the antiproton-to-proton ratio 
in Figure [7] which certainly tighten the constraints. 
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FIG. 6: The antiproton-to-proton ratio p/p for different 
branching fractions with m x — 300GeV and BF=110. 
Br(xx t + t~) + Br(xx ~ ► bb) = 1 in our model. In this 
figure, Br(xx ~ > t + t~) ~ 0.1 is typical for the usual MSSM. 
All the other curves, however, cannot be found in the usual 
MSSM but are possible in the uplifted Higgs region discussed 
in this paper. 



For a reference, we also showed in Fig. in addi- 
tion to MED model we have been using throughout the 
paper, the MIN and MAX models to indicate the fur- 
ther uncertainties arising from the propagation parame- 
ters [7l| . MAX models would have rather stringent 
constraints from PAMELA p/p data while the antipro- 
ton productions for MED and MIN parameters can still 
be within the observed bound. Despite those uncertain- 
ties, one can see that the our MSSM scenarios with the 
O(10)% bb channel can be compatible with the currently 
available p/p data while producing the positron excess. 



FIG. 7: The antiproton-to-proton ratio p/p including the 
background contributions for different propagation param- 
eters with m x = 300GeV, BF=110, Br( X x -> r + r~) = 0.9 
and Br(xx -> bb) = 0.1. 
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C. Gamma Ray Constraints 

The 7-ray data can also potentially exclude some of 
uplifted Higgs regions because our scenarios expect to 
give a significant 7-ray flux from r decays |63l - [65l l88l - [91| . 
The 7-rays, being chargeless, travel undeflected, and the 
flux from DM annihilation is given by 



d$ 7 



1 r Q p 2 Q{av) dN. 

2 Anml ~d~E- 



7 J(AQ)An 



with 



J(AQ) 



An 



J(tp)dQ 



(12) 



(13) 



FIG. 8: The gamma ray fluxes from the Galactic Center 
due to DM annihilation for the solid angle Afi = 10 sr. 
m x =300GeV, BF=110, BR(xx -»■ ^ + ^~)+BR(xx ->■ bb)=l. 
We also show, for reference, the HESS data along with its 
fitting function. NFW and Iso refer to NFW and isothermal 
halo profiles. 



and 



J(V0 



1 



r&P Q Jo 



p\r)dm 



(14) 



AS! 



where ip is the angle between the line of sight and the 
direction of the galactic center, r 2 = I 2 + pg — 2lp & cos ip 
and Afl = 2ir(l — cos ijj) is the observed region of the 



8 



DM: Br(xx->t*t")=0.9 ■ 
DM: BrfcOt-wV )=".7 
DM+BG: Br(xx^tV)=0.9 ■ 
DM+BG: Br(xx^rt>0.7 

} Galactic BG 
wExtra-galactic BG • 
J 




I (GeV) 



FIG. 9: The gamma ray fluxes for 10° < \b\ < 20°, 0° < / < 
360°. m x = 300GeV, BF=110, BR(xx ~> ^ + t")+BR(xX ->■ 
bb) = l. Data are from FERMI and the isothermal profile is 
used. The background (BG) that is added to the DM signals 
is the sum of Galactic and Extragalactic backgrounds, each 
of which is plotted here. 



sky. Stringent 7-ray constraints on our DM annihilation 
scenarios with the dominant tau channel arise from the 
prompt 7-rays due to the fragmentation of the annihila- 
tion products. 

The previous gamma ray data even before FERMI 
such as HESS [92[ already 



even 
offer tight 



constraints. 



Figure [8] illustrates the HESS data from the Galac- 
tic Center, its fitting function E 2 d$ 7 /dE = 2.3 x 
10" 12 (£ , /Tey)" 2 - 25 [r e y- 1 cTO- 2 s- 1 ], as well as a num- 
ber of predictions from DM annihilation. We have plot- 
ted the gamma ray flux from the Galactic Center (GC) 
for our canonical example due to the hadronic process 
calculated by DarkSUSY [93| where Pythia is im- 
plemented. For comparison, we also plot the 7-rays for 
a smaller tau channel branching ratio, which makes the 
spectra softer as expected. 

The DM predictions towards the GC within the re- 
solved direction can be strongly dependent on the halo 
density profiles. The Navarro Frenk White (NFW) pro- 
file [6a|, for instance, has a density that is strongly 
peaked towards the GC and leads to 7-ray fluxes in ex- 
cess of observations. The isothermal profile, which has a 
central core, predicts fewer 7-rays and is in better agree- 
ment with data. For this reason we have been using the 
isothermal profile throughout the paper. 

An unprecedented plethora of all-sky gamma ray data 
from FERMI offers further constraints on dark matter 
annihilation scenarios. 

The prompt gamma ray contributions along with the 



the preliminary FERMI data in the mid latitude 10 < 
b < 20 95ij extending the energy range published in Ref 
[j| are shown in Fig. [SJ This can give an upper bound on 
the gamma ray flux from the dark matter annihilations 
in our scenarios. 

To be more restrictive, we also show, for illustra- 
tion, the background contributions (which however have 
large uncertainties as well as contamination from point 
sources). For a simple estimation of the background, 
we added the Galactic background contributions (the 
sum of contributions from inverse Compton scattering, 
Bremsstrahlun g an d tt ) of the conventional GALPROP 
model in Ref. [96| and the isotropic diffuse extragalac- 
tic 7-ray background from FERMI parametrized by the 
power law <J>^ G oc E~ 2A1 [93] ■ We can see, despite the 
the uncertainties in the background estimation, that the 
dark matter annihilation scenarios in the MSSM with 
the electroweak scale mass can be still compatible with 
the current 7-ray data. 

It will be of great interest to see what insights into 
supersymmetric model building can be obtained from 
upcoming cosmic and 7-ray data with reduced system- 
atic errors and improved background estimation, par- 
ticularly when combined with upcoming data from the 
Large Hadron Collider. 



IV. CONCLUSION/DISCUSSION 

We have illustrated that dark matter annihilation sce- 
narios in the MSSM can be viable candidates to explain 
all current cosmic ray observations without necessarily 
introducing any additional degrees of freedom. In par- 
ticular, we have studied pseudoscalar Higgs s-channel 
resonance in the uplifted Higgs region to obtain boosted 
leptophilic annihilation cross sections which can explain 
PAMELA data while not in conflict with any other data 
sets. 

We here mention several directions for possible future 
improvement of our analysis. The positron fraction con- 
straints mainly come from the highest energy bins of the 
data with large error bars and it should improve by the 
forthcoming data from PAMELA and AMS-02 with bet- 
ter controls of the systematic/statistical errors. A proper 
treatment of the charge dependent solar modulation ef- 
fects (which we did not consider) could account for dis- 
crepancies in the low energy data (E < lOGeV) among 
different experiments. In addition to the positron frac- 
tion which also suffer from the electron background esti- 
mation uncertainties, the absolute flux of positrons such 
as those from AMS-02 would give more definite probe of 
the underlying physics. 

Our studies assumed the same value of the boost factor 
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for all species of cosmic rays. However, this is not neces- 
sarily the case. In particular, the antiprotons we observe 
originate from a large region of the halo in contrast to 
the positrons which are produced locally (say within a 
few kpc). Hence, in principle, the antiproton flux could 
be less boosted than the positron flux, for instance if the 
positrons are boosted partially by the local dumpiness 
(5J-[59|. Such relaxation of the antiproton constraints 
could be important for the annihilation scenarios within 
the framework of the MSSM, most of which are excluded 
due to antiproton overproduction. 

We presented our discussion primarily assuming an 
isothermal profile for the galactic halo. In fact a more 
cuspy profile towards the center, such as an NFW profile, 
was shown to be severely constrained (63M65| . Currently 
there is a great deal of uncertainty regarding the den- 
sity distribution in the inner regions of galaxies, though 
there exist observations of flattened cores [98| in spiral 
galaxies. More precise observations as well as detailed 
simulations including the proper treatment of the gas 
physics will be required for more realistic modeling of 
the halo profiles. 

Further studies of 7-ray constraints beyond the anal- 
ysis of this paper would be warranted. We refer the 
reader to, for instance, Ref. (63M651 |99| for gamma ray 
constraints from additional regions of the sky and halo 
profiles. We focused on gamma rays from the hadronic 
processes which are characteristic of DM annihilation 
scenarios in which r channel dominates; the consider- 
ation of additional effects such as the inverse Compton 
scatt ering inc luding those to the extragalactic gamma 
rays 100Ml02| could give additional constraints depend- 



ing on the parameter ranges of interest. 

More stringent constraints could however come from 
the particle physics rather than from the astrophysics 
once we have a concrete particle physics model. For 
instance, for the uplifted Higgs scenario, the tuning of 



the heavy Higgs decay width Ta/itia < 0(1O~ 3 ) implies 
y T < O(I0 _1 ), which needs to be checked with flavor 
physics constraints. We also could in principle consider 
a large in an analogous manner to a large y T by ad- 
justing the slepton, squark masses and the phase of the 
gluon mass [60(. If fj, final states could be significant by 
such (possibly fine-tuned) adjustments, then the 7-ray 
constraints could be relaxed relative to the r dominant 
scenarios considered in this paper [35l [36l lo3l l&ij even 
though here again the particle physics constraints (e.g. 
the flavor changing neutral currents) would need to be 
carefully checked. 

In the vast MSSM parameter space, it is worth search- 
ing for other regions previously missed that could explain 
all the cosmic ray data. The uplifted Higgs region it- 
self deserves further study. Even though our cosmic ray 
analysis is based on the properties of the uplifted Higgs 
regions, we kept our analysis fairly general so that a sim- 
ilar study could be applicable to any other parameter re- 
gions with similar properties. The unprecedented wealth 
of data expected from upcoming astrophysical and ter- 
restrial experiments (e.g. AMS-02 and LHC) could well 
provide further motivation towards a fuller exploration 
of the MSSM parameter space in the coming years. 
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